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Abstract 

The role of classical noise in quantum walks (QW) on integers is in- 
vestigated in the form of discrete dichotomic random variable affecting 
its reshuffling matrix parametrized as a SU2)/U(l) coset element. Anal- 
ysis in terms of quantum statistical moments and generating functions, 
derived by the completely positive trace preserving (CPTP) map govern- 
ing evolution, reveals a pronounced eventual transition in walk's diffusion 
mode, from a quantum ballistic regime with rate 0{t) to a classical diffu- 
sive regime with rate 0{y/t), when condition {strength of noise parame- 
ter) 2 x {number of steps) = 1, is satisfied. The role of classical randomness 
is studied showing that the randomized QW, when treated on the stochas- 
tic average level by means of an appropriate CPTP averaging map, turns 
out to be equivalent to a novel quantized classical walk without random- 
ness. This result emphasizes the dual role of quantization/randomization 
in the context of classical random walk. Keywords: Quantum walk, ran- 
domness, CP map, quantization 



1 Introduction 



Quantum Walk on Z : Recapitulation. The essential feature of a QW on a line 
is the promotion of the mathematical correspondence: left/right—^ head/tails, 
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between walker's move directions and coin's sides, to a dynamic interaction 
among two physical quantum systems PQ- [3]. This is realized by introducing state 
Hilbert spaces H w =sp&n(\m)) me z ~ C z and H c =span(|+), |— )) « C 2 , for the 
quantum walker and coin systems respectively. In H w the algebra Euclidean 
group 150(2) with generators {L, E + , £L} satisfying relations [L, E±] = ±E±, 
[E+,E—] = 0,is represented by the step operators E±\m) = \m ± 1} and the 
position operator L\m) = m\m). Introduce also the coin space projection 
operators P± = |±)(±|. One step of a CRW is described by the unitary 
V c \ = P+ ® E + + P- ® E_, acting on the total space H c ® H w , and via its 
conditional action on walker states, realizes the walker's move. Explicitly, coin 
p c and walker p w density matrices, initially in a product state p c ® p w , interact 
unitarily via the map p c ® p w — > V c \p c <8> PwV^ =AdV c \(p c <E) p w ) (written alter- 
natively in terms of the adjoint action AdX(p) := X(p)X^), and then de-couple 
by an unconditional quantum measurement of the coin subsystem, realized by 
the partial trace Tr c , i.e. V c \p c ® p w Ki ^•Tr c (T4ip c <£> PwV^). This constitutes 
a dynamic realization of the coin tossing process. The diagonal elements of 
the resulting walker's density matrix £ c i(p w ) =Tr c (AdV^i(p c ® Pw)), realizes the 
CRW on the integers [5],[B]. 

Quantization of CRW is conceived as the incorporation in H c of an addi- 
tional unitary operator U, the coin reshuffling matrix, in two different modes: 
either via the map AdV^i — >AdV^ioAdt/ €5 1, resulting into the U- quantization 
rule, or via the map AdT4i —>AdV c i o £ c :—AdV c \ o J^. A^AdL^ ® 1, for given 
ensemble {Aj, Ui} i of reshuffling matrices U\, forming a CPTP0], quantization 
map £ c '-=J2i K^dUi, resulting into the E- quantization rule [BJ. 



2 Noisy diagonal QW 

Consider the unitary 'walker-coin' operator V(s) = V c \U{s) ® l w , or explicitly 
V(s) = P+U(s)®E++P-U(s)®E-, (1) 

where the reshuffling matrix U(s), according to the Appendix, is taken to be 
an element of the coset SU (2) /U (1) » CP. Referring to (|2"2"j) . the coordinate of 
the complex projective variable 9 € CP, is chosen to be a dichotomous random 
variable (rv) 9 = se, e > 0, with values # = ±e, where s = ±1, with equal 
probabilities. The reshuffling matrix 

U(s)=Af(_ 1 se (2) 

with Af = ^^ +t 2 1 provides the unitaries V(s = ±1) that apply at random, with 
probability q± =1/2 (unbiased walk), so that the 'walker-coin' density operator 
evolves from the A-th to (N + l)-th step as 

PN+i = lYl V(s)p N V(s)l (3) 

1 s=±l 



2 



Rewrite this as p N+1 = £(p N ) = (V(s)p N V(s)*) = \ Y, s =± v q( s )PNV q (sY = 
^AdV(s)pN, where £ :Lin(i? c ® H w ) -^Lin(H c <g> H w ) stands for a CPTP map 
acting on the composite coin- walker system which implements the ensemble av- 
erage of the random unitary evolution. This density matrix iteration provides 
the sequence {pjv}jv=i! which subsequently provides, for any observable oper- 
ator A c eg) B w on H c <g> iJ w , a sequence of quantum moments {/ijv}jv=i! as 
Pn+i =Tr(pjy +1 A c ® B w ) = (A c <g> B W ) N . Utilizing the state-observable (trace 
inner product) duality, we derive the moments 

p N+1 = Tr(£ (p N )A c ® B w ) = Tr(p N r (A c ® B w )) = Tc(f 1 i 1, *(A c (8 B w )), (4) 

either evolving the state £, or the observables by the dual CP map £ * . The latter 
has as Kraus generators, the hermitean conjugates of those of £ . By means of 
the dual map £* , moment iterations for various choices of observables A c Cgi B w , 
are investigated next. 

Let p N = E Q , ££{+,-} J2k,\ez(PN)k, a .\,i3 I a) (P\ ® (Z| , be the density ma- 
trix describing the state of the "coin-walker" system, where (pN)k,a,\,p are ele- 
ments labelling walker positions k, I € Z, and coin states a, /? G {+, — }. Assume 

initially a factorized state = ®Po > Po = |0) (0| > Po =diag(cos 2 7, sin 2 7), 
for 7 € [0, 7r). The resulting matrix valued recurrence relation implies that pn 
is also diagonal. Setting (/9/v)kk =diag(aAr/c, /3zvfc)> for diagonal components 
obtained by expectation values of the projectors, 

a Nk = Tr( j0N P+ ® P k ) = (P+ <g> P k ) N , (5) 
/3jvfe = Tr(p N P_ ® P k ) = (P+ <g> P k ) N . (6) 

leads to recurrence relations 

a N+ i k = N 2 (ajvfc-i + e 2 /3ATfe-i) , (7) 
Av+ifc = A/" 2 (6 2 ttJV t + i +^jvfc+i) i (8) 

for N = 0, 1, 2, . . ., k £ Z, with initial conditions ctoo = cos 2 7, /?oo = sin 2 7, 
"ofc = Pak = for fc ^ 0. 



3 Recurrence relations for moments 

The recurrence relations © are difficult to solve exactly, except when e = or 
1. We consider instead the quantum moments of the distribution of probability 
over fc-values ('walker positions') after N iterations. Accordingly, let A$ = 

EZ-oo ^NkkP, and B^ = EZ-oo PNkk* ,for p = 0, 1, 2, . . .. Then A%> 

is the p-th moment of the distribution of probability over k values, with coin 

(p) 

spin 'up,' and B£ is the p-th moment of the distribution of probability over k 
values, with coin spin 'down.' The sum A$ + B)fi is the total p-th moment of 
the distribution of probability over k values. The operator expression of these 
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moments and their combinations sffl = A$ + B^' , and D^' — A$ — B$ , 
read 

A$ = Tr(P + ®LPp N ),B^=Tr(P_®LPp N ), (9) 
D%> = Tr((P + -P_)®LP j0N ), (10) 
= Tr((P + + P_)®LPp N )=Tr(l c ®LPp N )=:Tr(LP^ w) ). (11) 

The resulting recurrence relations for the moments are also difficult to solve 
in closed form for general N and p, except small values of p. In this way we find 
that 

A<g> =\(1 + r N cos(2y)) , < = \ (l - r N cos(2-y)) . (12) 
For the higher moments, we find 

D$ = (1 - r N )/(l - r) S$ = rc OS (2 7 )[l - r N ]/[l - r] , (13) 

and 

D {2) = 2rcos(2j)(l-r N )/{l-r) 2 - Nr N cos(2-f)(l + r)/(l-r) (14) 
Sff = N{l + r)/(l-r)-2r(l-r N )/(l-r) 2 . (15) 

The second total moment i.e. the expectation value of the square of the 
position operator L, can also be written as 

:( 2 ) _ rjTj.ff 2„(w)^ _ /f 2 
1 



S%> = Tr(LVD^(L 2 ) N (16) 



^(2Ne 2 - 1 + e 4 - (1 - e 2 ) N+1 /(l + e 2 )^ 1 ). (17) 

There are two regimes of interest for the second moment, corresponding to 
two value ranges of discrete time steps N relative to the strength of randomness 
coefficient 1/e 2 . The first regime is when 1/e 2 ^> N 3> 1, and then 

(1 - e 2 ) N+1 /{l + e 2 )^" 1 = 1 - 2Ne 2 + (2N 2 - l)e 4 + 0{{Ne 2 f), (18) 

so that 

4 2) = ^ 2 + ^((^ 2 ) 3 )- (19) 

This is the 'ballistic' or 'inertial' regime, where the growth rate of the second 
iment v 

(qw) m- 



moment with N is J (^L 2 ^ ~ O(N), that is characteristic of a quantum walk 



The second regime is when JV > 1/e > 1, and then 



S^ = ^{l + 0(l/(Ne 2 )}, (20) 



4 



i.e. y \L 2 j ~ showing the typical growth rate of a CRW or diffusion 

process. Note that, however small the e 2 , eventually N becomes so large that 
the second regime is reached. The changeover occurs for Ne 2 s» 1. Thus we see 
that for very small e 2 , the process behaves — as far as the second moment is 
concerned — initially like a QW, but eventually like a CRW. This is precisely 
the behavior observed in numerical simulations of a variety of quantum walks 
with unitary noise [7] and in previous studies e.g. [8]-[llj. of which the present 
process may be considered a simple, special case, with the advantage that it is 
amenable to a complete mathematical analysis. 



4 Randomization and Quantization of Walks 

Let us return to the evolution equation ([3]) for the total system density ma- 
trix. The following elaboration is possible; factoring the evolution operator as 
V(s) = V c \U(s) (8 l w makes explicit that the stochastic averaging amounts to 
averaging with respect to the matrices U(s = +), U(s = — ) with probabili- 
ties p± = |, which is implemented by a CPTP map £ c — \ S s=± Ad(i7(s) <X> 
l w ), acting non trivially only on H c . Explicitly we obtain pn+i = £(pn) = 
V cl £ c (p N )V^=AdV cl o£ c ( PN ). 

Noticeable, the averaging map £ =AdV c \ o £ c can be identified to be itself 
the CPTP evolution map of another QW. This new QW will have as its "walker 
space" the total Hilbert space of the initial walk i.e. H t := H c <Ei H w . In 
addition to the operator sum decomposition of the map £ defined following 
[3l it could also be described by a unitary dilation matrix Y in a extended 
space -ffaux ® H c ® H Wl where now an auxiliary vector space -ff aux should be 
attached to the initial space H t , as an additional "coin space". Such particular 
QW-like unitary dilation for £ is given in terms of the unitary operator Y = 
Y^,s=± P sQ ® v ( s )> 01 Y = Y ciQ ® It, acting on iJ aux ® H c ® H w . Under fairly 
general conditions, for CP maps with unitary Kraus generators (and £ provides 
an example), such a QW-like unitary dilation it is always possible (cf. [12], for 
a proof and related discussion). Explicitly for £ we have that 

p N+ i = £{pn) = Tr aux AdY(p aux (g> p N ) = Tr aux AdY cl o (AdQ (g) l t )(p aux ® Pn)- 

(21) 

Let p aux = I ft) (ft| be auxiliary system's state with \a) a basis vector. The 
non- unique unitary Q is chosen so that q± — (±| Q o Q* \a), i.e. 

^f^_ | . If o i = I matrix Q is the 7r/4 rotation matrix. 

In conclusion, the QW on integers with evolution map £ =Tr c AdV^Ad£/, 
randomized by un-bias dichotomous noise, when treated at the stochastic aver- 
age level with map £ =Ady c io£ c , becomes equivalent to a novel non-random QW 
with map £ =Tr aux AdF c iAdQ. The latter QW, via iterative relations obeyed by 
quantum moments of its CP evolution map, exhibits an eventual transition from 
a ballistic regime of fast spreading rate to diffusive regime of slow spreading 
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rate. Similar transition phenomenon occurs and quantization/randomization 
interplay is exhibited when randomization via continuous random variables is 
introduced. |13j 

Appendix: Lemma 1. The probabilities p, determining the CP map of QW 

£(Pw) = J2i=±Pi E iPw E i> are S iven b y Pi = (i\W o W* |c), where W is an 
element of the coset SU(2)/U(\) ~ CP coset corresponding to the reshuffling 
matrix U, and the basis vector |c) is determined by the initial state of the coin. 

Proof. The Hadamard or element- wise product (AoB)ij = AijBij of matrices 
A,B £ M{C N ), is invariant (DAo (DP.)* = AoB*) under diagonal D unimodu- 
lar (\Du\ — 1) matrices. Let V q = V C \U <E> It, be a unitary dilation of CP map £, 
viz. £{p w ) =Tr c (V q p c (g) p w Vj) = J2i=±Pi E iPwEj, where pi =Tr(P i J7p c ?7 t ) = 
(i| U o U* |c) , for Pc = |c) (c| . If U = DW be the SU(2)/U(1) coset decompo- 
sition of U, for D £ U(l), we obtain p l = (i\ W o (W)* |c) . If D = exp(i</x7 3 ), 
W (77) =exp(77(j + — ?7*cr_), we obtain 



1-M 2 ^ \ 1/1 



^f) = v / % = 1 9* 1 ' W 



where z = ^ = , e , € CP. In terms of p := (1 + lAr)" 1 = 1 - M , the 
bi-stochastic matrix W o (W)* = pl+(l — p)o\ reads 

w ° lw) " = 77w(m> '?)■ (23) 
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